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Abstract : We analyze equilibrium properties and adsorption desorption 
phase transition behaviour of a linear semiflexible copolymer chain under 
constrained geometrical situation on square lattice in a good solvent. One 
dimensional stair shaped line imposes geometrical constrain on the chain. 
Lattice model of fully directed self avoiding walk is used to model the chain, 
semiflexibility of the chain is accounted by introducing energy barrier for 
each bend of the chain. Exact expression of the partition function of the 
chain is obtained using generating function technique for the cases, viz. (i) 
constrained copolymer chain is in the bulk, (ii) constrained copolymer chain 
interacting with an impenetrable flat surface, (iii) constrained copolymer 
chain interacting with constraint itself and (iv) general expression of the 
partition function of the copolymer chain, interacting with a flat surface and 
geometrical constraint (stair shaped line). We have compared bulk proper- 
ties and adsorption desorption transition behaviour of a linear semiflexible 
homopolymer chain without constraint to the case when the chain is con- 
strained. 

PACS Nos.: 05.70.Fh, 68.18.Jk, 68.47.Pe, 36.20.-r 

1 Introduction: 

Biomolecules {DNA & proteins) live in the crowded, constrained regime. 
Such molecules are soft object and therefore can be easily squeezed into 
spaces that are much smaller than the natural size of the molecule in the 
bulk. For instance, actin fllaments in eukaryotic cell [1] or protein encap- 
sulated in Ecoli [2] is found in nature are examples of conflned biomolecules 
serves the basis for understanding numerous phenomenon observed in poly- 
mer technology, bio-technology and many other molecular processes occur- 
ring in the living cells. The conformational properties of single biomolecules 



1 



have attracted considerable attention in recent years due to developments in 
single molecule based observations and experiments [3, 4, 5, 6, 7, 8]. Under 
confined geometrical condition, excluded volume effects and effect of geomet- 
rical constraint compete with entropy of the molecule. Therefore, constraints 
can modify conformational properties and adsorption desorption transition 
behaviours of macromolcculcs. 

Behaviour of the flexible polymer molecule in a good solvent, confined 
to different geometries have been studied for past few years [9, 10]. For 
example, Brak and his coworkers [11] used directed self avoiding walk model 
to study behaviour of flexible polymer chain conflned between two parallel 
walls on square lattice. However, in present investigation we have considered 
a linear semiflexible copolymer chain on square lattice constrained by a one 
dimensional stair shaped surface. For two dimensional space, surface is a line 
and copolymer chain is constrained by the such surface. The chain is made 
of two type of monomers A and B distributed along the back bone of chain 
in a sequence A^B~A — B.... Therefore, it is a sequential or alternating 
copolymer chain. We have considered sequential copolymer chain because 
it serve as a paradigmatic model of actually disordered macromolecule e. g. 
proteins. 

To analyze the effect of constraint, we have used fully directed self avoid- 
ing walk model introduced by Privmann et. al [12] and used generating 
function technique to solve the model analytically. The results so obtained 
is used to discuss the behaviour of homopolymer chain in constrained geom- 
etry and also to compare results obtained for conformational properties of 
the polymer chain for the case, when chain is in the bulk, in absence of one 
dimensional stair to the case when there is stair to constrain the chain. If 
constraint is an attractive surface, it contributes an energy (< 0) for each 
step of walk making on the surface. This leads to an increased probability 
defined by a Boltzmann weight ui — exp(— es/ZceT") of making a step on the 
surface (e^ < or > 1, T is temperature and ks is the Boltzmann con- 
stant). The chain gets adsorbed on the surface at an appropriate value of u 
or eg- Therefore, transition between adsorbed and desorbed phases is marked 
by a critical value of adsorption energy or uic- The crossover exponent 
at the adsorption transition point is defined as, Ng ~ N"^, where N is the 
total number of monomers in the chain while A^^ is the number of monomers 
adsorbed onto the surface. A comparison is also made to discuss effect of 
constraint on adsorption transition point. 

The paper is organized as follows: In Sec. 2 lattice model of directed 
self avoiding walk is described for a linear semifiexible sequential copolymer 
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chain in good solvent. In sub-section 2.1, we have discussed behaviour of 
the copolymer chain in the bulk in presence of one dimensional stair shaped 
surface. The results obtained in constrained geometry case is compared with 
those obtained in absence of constraint for semiflexible homopolymer chain. 
Sub-section 2.2, is devoted to discuss adsorption of copolymer chain on a 
one dimensional flat surface in presence of constraint and results obtained 
is compared for the case when there is no constraint near the homopolymer 
and copolymer chains. In Sub-section 2.3, adsorption of copolymer chain on 
constraint is discussed. While in sub-section 2.4, general expression of the 
partition function of copolymer chain which is interacting with constraint 
and fiat surface is discussed. Finally, in Sec. 3 we summarize the results 
obtained. 

2 Model and method 

A lattice model of fully directed self-avoiding walk [12] is used to model a 
semiflexible alternating copolymer chain in two dimensions and generating 
function technique has been used to calculate partition function of the copoly- 
mer chain in presence of an impenetrable surface (line) having shape like a 
stair (as shown schematically in figure 1). Since copolymer chain is fully 
directed therefore, walker is allowed to take steps along +x, +y directions 
on the square lattice. The stiffness of the chain is accounted by introduc- 
ing an energy barrier for each bend in the walk of the chain. The stiffness 
weight k = exp{—f3eb) where (3 = {kBT)~^ is inverse of the temperature and 
€(,(> 0) is the energy associated with each bend of the walk of copolymer 
chain. For A; = 1 or e;, = the chain is said to be fiexible and for < /c < 1 
or < €6 < oo the chain is semiflexible. However, when — >■ oo or /c — >■ 0, 
the copolymer chain is like a rigid rod. 

Directed self avoiding walk model can be solved analytically and therefore 
gives exact value of conformational properties and adsorption transition point 
of the chain. However, directed walk model is restrictive in the sense that 
bending angle has same value i. e. 90° (for present study) for each bend and 
directedness of the walk introduces certain extent of stiffness in the polymer 
chain because all directions of the space are not treated equally by the walker. 

The partition function of a semiflexible sequential copolymer chain made 
of two type of monomers (A & B) can be written as. 
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Figure 1: In this figure polymer chain is shown (I) in absence of constraint 
(one dimensional stair shaped surface) while in (II) there is constraint. Walk 
of the chain starts from point O and directed along both the directions of 
space. 



Z{k,x„x,) = j:l- E x,^/'x,^V^ (1) 

all walks of N steps 

where, Nj, is the total number of bends in a self avoiding walk of steps 
(monomers), Xi and X2 is the step fugacity of each of the A and B type 
monomers, respectively. For the sake of mathematical simplicity we assume 
here and onwards Xi = X2 — x. Method of analysis discussed in this paper 
can also be extended to the case when xi ^ x^- 

We consider an impenetrable line, located in the x — y plane (as shown 
schematically in the figure 1) has shape like a stair, one end of the copoly- 
mer chain is grafted to origin O. Since stair is impenetrable, therefore, it 
constrains the polymer chain. We are interested in analyzing effect of this 
constrain on the conformation properties of the copolymer chain in the bulk 
and also in calculating adsorption desorption phase transition point of the 
copolymer chain in presence of constraint. 

2.1 Copolymer chain in the hulk 

The partition function of a sequential semifiexible copolymer chain in the 
bulk is calculated using the method of analysis discussed by Mishra et al, 



4 



[13] and components of the partition function Z[k, x) of the copolymer chain 
in the bulk is written as follows: 

Xo{k,x) ^x + x'^k{l + Y + kXo) (2) 

where, Xn{k,x) (here and onwards n=0,l,2,3,4...) is the sum of Boltzmann 
weight of all walks with first step along +x-direction and n step distant from 
the stair while Y(k,x) is sum of Boltzmann weight of all the walks of the 
copolymer chain having first step along +y direction. Above equation can 
be re-written as, 

Similarly, we have 

"V // \ 2 n n'^^^ kx^ . i 2 in n ^ Xa/" 

y\fi\k,^X) = X-\-X -\-..-\-X -\-X — oT"o \~(rvX-\-kX -\-...-\-kX -\-X - oi o)^ 

1 — x^A;^ 1 — x^A;^ 

(4) 

and substituting n — oo, we recover X^{k, x) — X{k, x) [13], 

Xoo{k, x) = + 5 

1 — X 1 — x 

Where, Y' is the other component of the partition function in absence of con- 
straint. But, another component of the partition function of the semiflexible 
copolymer chain in presence of constraint is, 

Y{k, x) ^ X + x{x + kXi) + x'^{x + kX2) + x'^{x + kX^) + (6) 

so that, 

x{l — x'^){l — x'^k'^) + kx'^il — k'^x'^) + kx'^il — x){x + kx"^) 
" (1 - x){l - x^){l - k^x^) - kH\l - x^k^) - k^x\l - x) ^ ^ 

Thus, the partition function of the copolymer chain in the bulk is written as, 

Z{k,x) = Xo{k,x) + Y{k,x) (8) 

such that 

_ x{2 + (-1 + 2k)x + (-2 + k- 2k'^)x'^ + (l-2k + 2k'^ - 2k^)x^) 
^ ' ~ 1 - X - (1 + 2A;2)x2 + (1 + A;2)x3 + k^x^ ^ ' 
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From singularity of the partition function Z{k,x), we obtain critical value 
of step fugacity, Xc, needed for polymerization of an infinitely long linear 
semiflexible copolymer chain. Since, we have considered Xi — X2 — x, we 
can treat this result for semiflexible homopolymcr chain i. e. both the 
monomers of copolymer chain are identical. For, A; = 1, we have a;c(with 
constraint)=0. 532089 while ,Xc(without constraint)=0.5 [13] for flexible ho- 
mopolymcr chain for fully directed self avoiding walk model on square lattice. 
Variation of Xc with k in presence as well as absence of stair (constraint) is 
shown in figure 2. 

Persistent length, Ip, of the chain can be calculated using definition of 

Mishra et. al [13], lp{= ,dLogz(k,x) )■> ^ average distance of the polymer chain 

^ dk 

between its two successive bends. The variation of the persistent length with 
k is shown in figTirc 2. A comparison of Ip is also made for the case when 
there no constraint [13]. 

2.2 Adsorption on one dimensional flat line 

In order to analyze the effect of presence of constraint on adsorption transi- 
tion point, we consider an impenetrable line located at x=0 and calculated 
components of the partition function of the sequentional copolymer chain, 
which is interacting with impenetrable flat line (while there is no attraction 
between monomer and constraint) using method of analysis used by Mishra 
et. al, we have component of the partition function along the surface as. 



S^{k, uji, 002, x) = si(l + S2 + kXi) + 8182(1 + 82 + kX2) + SiS2(l + S2 + kX^) + . . . (10) 

here A type monomer of the copolymer chain is grafted at O, we have Si(= 
(jJix), S2 = {u!2x) is the Boltzmann weight of A and B type monomers of 
the copolymer chain with flat surface, respectively. While Xn{k,x) is the 
component of the partition function perpendicular to the flat surface i. e. 
along The value of Xn{k,x) wc have used from Eq. (4) to calculate 

partition function, Z'l'{k,U!i,(x!2,x), of the copolymer chain (provided 8182 < 
1) as follows, 

S^(k u 00 x) - Sikx{l + 82) ^ SiS2kx'^{l + Si) 

1 - S1S2 {1 - Si82){l - 8i82X^) {1 - 8i82){l - 8i82X^) 

where. 
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Figure 2: This figure compares the value of conformational properties of 
semiflexible homopolymer chain when there is constraint (an impenetrable 
stair shaped surface) to that of the case if there is no constraint. We have 
also compared the value of adsorption desorption transition point of the 
semiflexible homopolymer chain on a flat surface in presence and absence of 
the constraint. 



sikx{l+S2x){x+kx'^) 

and 

/ sik'^x{l+S2) I siS2k'^x'^{l+si) , si{l+S2x)k'^x^ \ 
^ {l—SlS2){l—SlS2X'^) (1— SlS2)(l~"SlS2X^) {l—siS2X'^){l—k''^x'^) ^ 
Therefore, partition function of the copolymer chain with A type monomer 

grafted at O can be calculated as, 



Z^{k,u!i,uj2,x) = S^{k, 001,002, x) + Xo{k,x) (siS2 < 1) (12) 

Z^{k, a;,, a;., x) = ^^^^^ _ ^ _ + ^f^2^^2 + + ^.)^3 + ^4^4))) (^i^^ < 1) (13) 
where, 

U3 = (-1 + SiS2){-l + SiS2x'^) 

We are interested in the singularity of the partition function. Therefore, 
comparing the denominator of the partition function Z^{k,oui,ou2,x) equal 
to zero, we obtain critical value of monomer surface attraction required for 
adsorption of copolymer chain, Ud = ^ on the fiat surface. Effect of 
constraint on adsorption transition of copolymer chain on fiat surface is an- 
alyzed through figure (3). Replacing si with S2 in Eq. (13) we could obtain 
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partition function for copolymer having first monomer of B type. The value 
of adsorption transition point remain unchanged when grafted monomer of 
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Figure 3: This figure compares the value of Ud required for adsorption of 
copolymer chain on the fiat surface for ujc2 — 0.5, 0.9 in presence and absence 
of the constraint. 



2.3 Adsorption of copolymer chain on constraint (one 
dimensional stair) 

The stair shaped one dimensional line (surface) is located in the x — y plane 
as shown schematically in figure (1). The copolymer chain is interacting 
with this surface. While, fiat surface (line) located at a; = is neutral 
for copolymer chain. The components of the partition function of surface 
interacting copolymer chain can be calculated using method discussed above, 
as. 



SoA{k,UJ3,UJ4,x) = Ss + S3S4k{l + YA + S3k)+sls4k'^{l + YA + S4k) + slslk^{l + YA + Ssk) + . . . 

(14) 
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where, S3(= oo^x), S4(= U4X) is the Bohzmann weight of interaction 
energy of A and B type monomers with stair shaped surface respectively. 
SnAik,uj2,uj4,x) [A indicates that first monomer grafted at O of the copoly- 
mer chain is of A type] is the sum of Boltzmann weight of all walks hav- 
ing first step along +a:-direction and n step distant from the stair while 
YA{k,U3,U4,x) is sum of Boltzmann weight of all the walks of the surface 
interacting copolymer chain having first step along +y direction. In the same 
fashion one can define Snsik, 003, 004, x) and YB{k, 003, 004, x). We assume that 
first monomer of copolymer chain is of A type therefore for S3S4 < 1, we 
have, 

S,A{k,U,,U4,x) = Y^J^ + (15) 

also, for odd n—p, 



and for even n—q, 



(16) 



The component of the partition function of surface interacting semiflexible 
copolymer chain along y-axis is, 

YA{k, UJ3, UJ4, x)^x + x{x + kSiB) +x'^{x + kS2A) + x^{x + kS^s) + ■ ■ ■ (18) 

substituting the value of SqA{k, uj3,uj4,x) and SpB{k, us,U4,x) in Eq. (18) we 
get, 

a;(l - x^)(l - S3S4k^) + kx'^{l - ^^5354) + kx^{l - x){s3 + ^5354) 

Ya[IC,UJ3,UJ4,X) - ^^^^ _ ^2^^^ _ _ ^2^2(1 _ ^35^^2) _ k^g^s^X-'f^l - x) ^ ' 

Thus, we calculated partition function of copolymer chain interacting with a 
stair shaped impenetrable surface as, 

Z''{k,UJ3,U4,x) = SoA{k,UJ3,UJ4,x) + YA{k, LJ3, LJ4, x) (20) 

therefore, partition function of the copolymer chain which is interacting with 
an impenetrable stair shaped line is, 

7./, X N''{k,U3,U4,x) , . 

Z [k, UJ3, UJ4, x) = — — (21) 

D'^{k,UJ3,LU4:X) 
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where, 

N'^=x + kx^ —x^ — S3{2k^S4X^ — ( — l + x)^(l + x) + k^x{s4 + x — 2s4x) + 
k{—l + x){s4 + S4X + x^)) 
and, 

= k^s^s^x^ + (-1 + xfil + x) + k\.S3S4i-l +x)- x^) 
From singularity of the partition function Z^{k,ujo^,uj,i.x), we found that 
= k^]^^x^(-^+x+k^x^) °^ comparing D'Xk^U'i^Ui^x) = 0) required for 

adsorption of semifiexible copolymer chain and it is independent of bending 
energy or stiffness k of the polymer chain. This we have discussed through 
numerical results shown in figure (4) for ujc2 = 0.5,0.9&1.0. This is because, 
the polymer chain is directed towards the surface (constraint). Replacing cu^ 
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Figure 4: This figure is used to show that copolymer chain adsorption on 
constraint appears independent of stiffness of the copolymer chain. However, 
in the very stiff chain limit adsorption of the chain on contsraint may not be 
possible. 



2.4 General expression of the partition function 

We now derive expression of the partition function of the sequential semifiex- 
ible linear copolymer chain for the case when copolymer chain is interacting 
with fiat surface and constraint. The fiat impenetrable surface (line) located 
at a; = and another surface is constraint (an impenetrable stair shaped line 
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located along the line x = y; schematically shown in the figure (1)). We as- 
sume that A type monomer of copolymer chain is first monomer of the chain 
and one end of the chain is grafted at (0,0) or O. Let -S"^ is the component 
of the partition function having first step on to the flat surface, then, 

SA{k, (jJi,(jJ2, (jJ4, x) = Si + SiS2k{l + SiB + 52^) + S^S2A;^(1 + S2A + 52^) + • ■ ■ (22) 

Substituting value of SgA, SpB from Eqs. (16& 17) and using following 
relation we have general expression of the partition function, 

Z^^^{k, oji, 0J2, 003, 004, x) = SAik, oji, 002, 003, 004, x) + SQA{k, ^3, UJ4, x) (23) 
Z*^-''(A;,a;i,a;2,a;3,a;4,x) = ^i^^—t^ +1*4 + ^5 - — + — (siS2 < 1; < 1) (24) 

1 - S1S2 U7 Ug 

where, 

ks^{l+ks4)six{l+S2x) 
''4"(-m^S3S4)(-l+siS2X^) 

ksi{l+S2)x{l+S2x) 

""'^ {-1+SIS2){-1+SIS2X^) 

Uq = S3(-{k'^(-l + S4)x'^) + k^S4x'^ - (-1 + xy(l + x) + kS4(-l + x^)) 

U7 = A;^S3S4X^ + (-1 + x)^(l +x)+ k'^{s3S4{-l +x) - x^) 

Us = {k'^six'^{—l — /c^S3S4( — 1 + x) + k^S3S4X + k{ — l + S3( — 1 + x))x + 

x'^){ — l — S2{l+X + Six)+SsS4{{ — l + SiS2){l + S2x) + k'^{l + S2 + S2X + SiS2x)))) 

and 

U9 = ((-1 + /C^S3S4)(-1 + SiS2)(-l + SiS2X^){k^S3S4x'^ + (-1 + xf{l + 

x) + A;2(s3S4(-i + x)- x^))) 

We call this relation (Eq. 24) as general expression of the partition func- 
tion because from this expression we can recover expression of partition func- 
tion discussed in above sub-sections by suitable substitutions. For example, 
replacing si, S2, S3 and S4 by x we have partition function of the copolymer 
chain in the bulk i. e. Z{k,x). 

3 Result and discussion 

We have considered a linear semiflexible sequential copolymer chain made 
of two type of monomers and fully directed self avoiding walk model is used 
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to model the copolymer chain on square lattice. Technique of generating 
function is used to calculate partition function of the chain in presence of 
constraint (an impenetrable one dimensional line) analytically. The variation 
of critical value of step fugacity required for polymerization of an infinitely 
long linear copolymer chain and persistent length under good solvent have 
been discussed. The results obtained with constraint is compared to case 
when there is no constraint. We have also compared the value of monomer 
surface attraction required for adsorption of semiflexible copolymer chain on 
a flat surface and constraint. The surface is a one dimensional line and walks 
of the chain is directed along both the directions of the space. 

We have found that adsorption of the homopolymcr chain on constraint 
occurs at a very small value of monomer and surface attraction. This is 
because, walks of the polymer chain is directed along the stair shaped surface, 
walks of the chain follows curvature of the constraint and major contribution 
in the partition function of copolymer chain is due to walks lying on the 
constraint. In the case of stair shaped impenetrable surface, critical value 
of monomer surface attraction required for the copolymer chain adsorption 
appears independent of bending energy of the chain (excluding the case of 
very stiff chains). While adsorption of copolymer and homopolymer chain 
on the flat surface is function of scmiflexibility of the chain. The stiffer chain 
adsorption occurs at a smaller value of monomer surface attraction than the 
flexible polymer chain. These features of adsorption on an impenetrable 
flat and curved surfaces are similar to the semiflexible polymer adsorption 
[13, 14, 15, 16, 17]. 

References 

[1] S Koster, D Steinhauser and T Pfohl 2007 J. Phys.: Cond. Matt. 17 S4091. 

[2] G Morrison and D Thirumalai 2005 J. Chem. Phys. 122 194907. 

[3] Structure and dynamics of confined polymers 2002 (ed. J J Kasianowicz et. 
al, Kluwer Academic Publishesr, Dordrecht). 

[4] S E Henrichson, M Misakian, B. Bobertson and J J Kasianowicz 2000 Phys. 
Rev. Lett. 85 3057. 

[5] D Nykypanchuk, H H Strey and D A Hoagland 2005 Macromolecules 38 145. 
[6] M Ichikawa, Y Matsuzawa and K Yoshikawa 2005 J. Phys. Soc. J. 74 1958. 



12 



[7] C H Reccius, J T Mannion, J D Cross and H G Craig-head 2005 Phys. Rev. 
Lett. 95 268101. 

[8] W Reisner, K J Morton, R Riehn, Y M Wang, Z Yu M Rosen, J C Sturm, S 
Y Chou, E Frey and R H Austin 2005 Phys. Rev. Lett. 94 196101. 

[9] P G de Gennes 1979 Scaling concepts in polymer physics (Cornell University 
Press, Ithaka). 

[10] E Eisenriegler 1993 Polymer near surfaces (World Scientific, Singapore). 

[11] R Brak, A L Owcarek, A Rechnitzer and S G Whittington 2005 J. Phys. A 
38 4309. 

[12] V Privman, G Forgacs and H L Frisch 1988 Phys. Rev. B 37 9897; Privman V 
and Svrakic N M 1989 Directed Models of Polymers, Interfaces, and Clusters: 
Scaling and Finite-Size Properties (Springer, Berlin). 

[13] P K Mishra, S Kumar and Y Singh 2003 Physica A 323 453. 

[14] P K Mishra, 2010 J. Phys.: Cond. Matt. 22 155103. 

[15] P K Mishra, 2010 Phase Tran. 83 47. 

[16] P K Mishra 2010 New York Sci. J. 3(1) 32. 

[17] P K Mishra 2009 Acad. Arena 1(6) 1. 



13 



